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Abstract. - We study theoretically the dimensional crossover from a three-dimensional elon- 
gated condensate to a one-dimensional condensate as the transverse degrees of freedom get 
frozen by tight confinement, in the limit of small density fluctuations, i.e. for a strongly degen- 
erate gas. We compute analytically the radially integrated density profile at low temperatures 
using a local density approximation, and study the behavior of phase fluctuations with the 
transverse confinement. Previous studies of phase fluctuations in trapped gases have either 
focused on the 3D elongated regimes or on the ID regime. The present approach recovers 
these previous results and is able to interpolate between them. We show in particular that 
in this strongly degenerate limit the shape of the spatial correlation function is insensitive to 
the transverse regime of confinement, pointing out to an almost universal behavior of phase 
fluctuations in elongated traps. 



In the recent years, one-dimensional (ID) ultracold atomic gases have been produced 
in very elongated magnetic [1,2] or optical traps [3,4], with such tight confinement in two 
transverse directions that the atomic motion "freezes" to radial zero-point oscillations. The 
equilibrium phase diagram of these trapped ID gases shows a rich behavior [5]: at low densi- 
ties the cloud behaves as a gas of impenetrable bosons ("Tonks-Girardeau gas" [6]), and for 
higher densities (corresponding to most current experimental setups), the cloud is a "quasi- 
condensate [5,7], characterized by suppressed density fluctuations and the same kind of local 
correlations as in a true condensate, but also by the lack of long-range phase coherence due to 
significant phase fluctuations. The latter die out continuously around a characteristic temper- 
ature T^, below which the quasicondensate turns into a true condensate. It was realized in [8] 
that such quasicondensates could also exist in very elongated, but three-dimensional (3D) 
traps: although atomic motion is possible in all directions, the lowest-lying excitations of 
these systems, which dominate the long range decay of the phase correlations, are ID in char- 
acter [9]. These 3D quasicondensates were observed in equilibrium [10] and non-equilibrium 
samples [11], and their coherence properties studied through Bragg spectroscopy [12,13] or 
matter- wave interferometry [14]. 
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Currently, although the limiting 3D or ID cases are well understood theoretically, the 
kinematic crossover from the "3D elongated" to the truly ID regime of transverse confinement 
have been less studied, despite being relevant to several recent experiments with "conventional 
traps" [1, 12], or to ongoing research on the manipulation of coherent atomic ensembles on 
micro-fabricated substrates [15, 16]. For the case of an infinitely long waveguide with uniform 
axial density, a many-body theory has been developed in [17] to describe this crossover (see 
also [18]). The influence of an axial trapping potential on the density profile and collective 
excitations was addressed in [18] by numerical integration of the transverse dynamics, in [19] 
through a separability ansatz, which is correct in the ID limit but not in the 3D case, and 
in [20], where solvable hydrodynamic models that reproduce the 3D and ID TF limits are 
introduced. 

In the present paper, we follow a different route and introduce an accurate analytical 
approximation for the "equation of state" of the quasi-lD gas, that connects smoothly to the 
3D to the ID mean- field regimes. In the local density approximation (LDA), we are then able 
to work out analytically the ID (integrated over radial degrees of freedom) density profile 
ni(z), from which many properties of the trapped quasi-lD gas can be calculated. As an 
appUcation, we investigate the influence of the transverse conflnement on phase-fluctuations 
in a quasi-lD geometry. Ref. [21] considers the related problem at zero temperature case in a 
box geometry, and investigate the 3D to ID crossover as the box cross-section is reduced. In 
the trapped case considered in this paper, quantum fluctuations of the phase are small away 
from the Tonks regime [5], so that we only consider thermal fluctuations. We calculate the 
spatial correlation function of the radially averaged atomic field operator ^, 



C^^\s) = Jdz{i>\z + s/2)^{z - s/2)) 



(1) 



at finite (but low) temperatures. This important quantity gives a global measure of phase 
coherence across the atomic cloud, and is also the Fourier transform of the momentum distri- 
bution (sec for instance [22]). The local density approach in this context is nothing else than a 
slowly- varying-envelopc approximation, applied to the long- wavelength excitations responsible 
for the fluctuations of the phase. We show that this turns out to be a good approximation even 
if T > T^, provided density fluctuations are small. More importantly, we find the remarkable 
property that the shape of the spatial correlation function (or equivalently the momentum 
distribution) of the quasi-condcnsatc is almost insensitive to the precise regime! of transverse 
confinement, pointing out to the universal character of phase fluctuations in very elongated 
traps. 

We follow the general method introduced in Ref. [5], assuming a weakly interacting gas well 
below the degeneracy temperature (see [23] and related discussion at the end of the paper). 
Then, density fluctuations are small and the equilibrium 3D density profile no in the trapping 
potential Vtrap(r) = muP^p^ /2 + V{z) is given by the solution of the usual Gross-Pitaevskii 
equation, even if the cloud is a quasi-condensate. Without loss of generality, we introduce the 
ID density ni{z) = J p no{p, z) obtained by integration over the transverse plane, and the 
radial mode /j^ through no(p, z) —\ f±{p,z) [^ ni{z). Assuming a sufficiently shallow axial 
conflnement and neglecting the density derivatives of ni , f± with respect to z obtains 



2M /i 



Mwip^ + Uni{z) 1 ]2= w.e.[ni(^)]. (2) 



Here, we have used cylindrical coordinates, denoting the transverse radius as p and the lon- 
gitudinal coordinate as z. The 3D coupling constant U is related to the s-wave scattering 
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length a through U — Anti^a/M, where M is the atomic mass. The local equilibrium chemical 
potential /ij c. depends on z through 

^il.c.[nl{z)] +V{z) ^ ^1, (3) 

with fj, the global chemical potential of the cloud in the trap. The relation between ni and 
fii.e., that includes the effect of transverse confinement, can be seen as an effective "equation 
of state" for the ID gas [18]. 

To obtain this equation of state, it is sufficient to solve Q locally, i.e. for each value of 
ni{z), or equivalently by considering a geometry with radial harmonic trapping, but homoge- 
neous Id density ni ~ N/2L, where 2L is the axial length of the system ("cylinder model", 
see inset in Fig. ^). The equation of state in this model was found in [18] by numerical inte- 
gration. Here we follow an alternative route, and derive the equilibrium properties by taking 
for /_L[ni] a Gaussian trial wavefunction, whose width k7^[71i] is a variational parameter, and 
by minimizing the chemical potential fJ-i.c.. This simple calculation yields the optimized width 
'*i'-L['^i] = +4:aniy^^, where a± — y/h/Muj± is the radial oscillator length, and the local 
chemical potential 

Mi.c.['^i] = hujj_\/l + Aarii. (4) 

The very good agreement of Eq. Q with the available numerical results [18], illustrated in 
Fig. nji, has been pointed out in [20] on a phenomenological basis. Here we show that this 
expression follows from the condition that /ii.o. is the lowest eigenvalue of the Gross-Pitaevskii 
equation. Note that this approximation does not correspond to a separability ansatz, as ni is 
in general a function of the axial coordinate z. In this respect, the method used here differs 
from the work reported in [17]. 

We now reintroduce the axial trapping potential, assuming for definiteness a harmonic 
form, V{z) — iAfw^z^. From the local equilibrium condition ^ and the equation of state 
(0}, one finds the density profile of the trapped gas. 



nr(.) ^ ^m-m 
16a nuj I 



ViL)-Viz) ^ ^ 
fku I 



£(l-z2)[«(l-52) + 4]. (5) 



Here we have introduced the parameter a = 2{ii/hLUj_ — 1) and z = z/L. The condensate 
length L is defined by the relation ^i,c.[ni{L) = 0] + ^AIoj^L'^ = /i, and is given explicitly by 

i=— V^, (6) 

where = ^/KJMlJ^ is the axial oscillator length. Using ni{z)dz — N, and the density 
profile lO, we obtain the equation for the key quantity a, 

a^a + 5f = (15xf. (7) 

The only parameter of the calculation, x — Naa±/a1, roughly gives the ratio of the interaction 
energy to the radial zero-point energy [18]. Numerical solution of Eq. Q is straightforward, 
and obtains the static properties of the condensate at any confinement strength. In the limit 
X ^ 5, the mean- field interaction dominate over the transverse confinement, and one recovers 
the well-known 3D TF result, a w asD = (15x)^^^. Conversely, if x ^ 5, the transverse 
motion is almost frozen and one finds a ~ aiD = (3x)^^^- The crossover between the two 
regimes occurs approximately for aiD = O-sd, giving a cross-over value Xcross = 5'^/^/3 « 3.73. 
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Our analytical results are very accurate even in the crossover region, as shown in Fig. 
where we compare Eq. ((SJ to a direct numerical solution of the GP equation. 

Equations (|5I6I7|I are the key results of this paper. The calculation of the density profile 
for any strength of the confinement in the transverse direction allows to deduce a number 
of interesting quantities. As an example, we focus in the remainder of the paper on phase 
fluctuations, a key phenomenon to understand the physics of ID gases. In general, phase 
fluctuations originate from the very large population of the ID excited states of the system, 
i.e. those in the low-energy energy range hujz < e ^ min{/i, fku^} (referred to in the following 
as the "axial branch" of excitations). Keeping with the local density approach, we analyze 
first these excitations in the cylinder geometry. The axial branch corresponds to excitations 
characterized by an axial wavenumber k, and no radial nodes [9]. Well below the degeneracy 
temperature [23] , the elementary excitations of the phase-fluctuating ensemble obey the same 
Bogoliubov-De Gennes equations than in the usual coherent ensemble [5]. We introduce 
the operators (f) and Sfi describing fluctuations of the phase and of the density, and their 
expansion on the set of axial plane waves, 5fi{z) — 5nkAk{p)e.^^^bk/ yf^jUiL -t- h.c, and 
(j)[z) = —i <f)ke^'''^bk/ yfiniL -\- h.c. The function Ak{p) describes the radial dependance 
of density fluctuations, and the operators Sfc, 6^, destroy and create one quasi-particle with 
wavevector k. The Fourier components of (j) and 6h read 

hk^ 

hwfSukAkip) = — no(p)0fc, (8) 



^k '/'/cno(p) = ( + Unoip) I SukAkip) - ^]^^J 



"o(/o)V_ 



Akjp) 
no{p) 



(9) 



For the axial branch of interest, the transverse envelope changes continuously from a flat 
profile in the 3D regime to the radial ground state in the ID regime [9]. In any case, the 
quantity Ak/no is almost flat; The corresponding spatial derivatives in (j^ are strictly zero 
in the ID limit, and of order {huj±/p)'^ <C 1 compared to the first term of the right hand 
side of lO in the 3D limit. We thus neglect their contribution, and average over transverse 
degrees of freedom to get rid of the remaining radial dependance [9] . The equations obtained 
in this way are of the usual Bogoliubov form, with an excitation spectrum given by spectrum 
ujf = ^{hk'^/2MY + clj^k"^, and the amplitudes Suk = ni{uJk/Lof Y^^ and 0^ = {ujf /AcjkY^^ . 
The longitudinal speed of sound is cl^{k) = UnoAk/M, and may depend on k through 
noAk = J p noAk for uj^ ^ . This has been suggested in [24] as a possible mechanism 
contributing to the decrease of the critical velocity in 3D elongated gases. The phase coherence 
properties are however not affected, since they are determined the phonon-like, lowest-energy 
modes with energy < {T /T^Y/'^hw^ ^ ^± [5,8], for which Ak is fc-independent regardless of 
the transverse regime [9]. 

We now include the effect of the trapping potential by introducing a local density profile 
ni{z) according to and ((SJ- In the above expressions, one than has to replace everywhere 
the chemical potential and density by their local values (see [13] for further details). Using 
(b^Sfe) ~ k^T/hhj^, the variance of phase fluctuations then reads 

A<A^(^,.) = mz) - ^t>{z')f) « f = (10) 

J0ni(z) L ni(z) 

with the relative distance s = z — z' and the mean coordinate z = (z + z')/2. The phase tem- 
perature in l|10|l is k-QT^ = h^ni{Q)/ML cx N{hijJz)'^ / p [5,8], and the phase coherence length 
is = LT^/T = h^ni{0)/MkBT. Together with the density envelope lO, the expression l{lT)|l 
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is sufficient to find the long wavelength behavior of the spatial correlation function [5] , 

C«(|) - ^ r dz V«i(z + s/2)ni(z-s/2)exp(-iA02(^,s)). (11) 

J — L 

This expression differs slightly from the one used in [13] in the treatment of the overlap term, 
defined here as y^ni(z + s/2)ni(z — s/2). The way we write it here yields better agreement 
with the numerical calculation of the correlation function in the 3D case (see FiglJb). 

Several comments can be made on this expression. First, already for T = AT^, the LDA 
expression (jllll agrees well with the numerical calculation of the correlation function based 
on Refs. [5,8]. This is shown in Fig|2K for the 3D case, and Fig|2l3 for the ID case. Second, 
the expression l|ll() . which depends on the dimensionless space variable s/L, is a universal 
function completely determined by two dimensionless parameters, Xj which controls the regime 
of transverse confinement and the functional form of the density profile, and T /T^^ which 
controls the magnitude of phase fluctuations. The third and most important conclusion is that 
the resulting correlation function is insensitive to a large extent to the transverse regime of 
confinement (in other words, to the value of x)- This is illustrated in Fig|2t), where we plot the 
correlation function for T = 4T^ and x — 100 (dotted), x = 1 (solid) and x = 10~^ (dashed). 
Despite the dissimilar transverse profiles, which correspond to very different experimental 
systems (see Table IJ), the axial correlation function are almost identical, showing almost 
exponential decay on a 1/e length scale ~ 1.54^0. This very weak dependance on x points 
out to the almost universal nature of thermal phase fluctuations in ultracold, very elongated 
trapped gases. Modifications of the functional profile rji(z)/ni(0), exphcitly present in Eq. 
Hll|l . are not significant. Rather, the effects of transverse confinement are almost entirely 
contained in the scaling variable oc ni(0). 

The approximate scaling identified here relies on the zero-temperature equation of state, 
Eq. 10}. Its validity requires that (i) axial density fluctuations are negligible, and (ii) that when 
going to the 3D regime, the presence of a normal cloud, mostly composed of 3D excited states 
with energy ^ hio^^ do not modify signiflcantly the density profile of the quasicondensate. 
For (i) to be true, one requires that T ^ Td, where Td = Nhujz is the ID degeneracy 
temperature [5,8]. Note this is always the case in the 3D case since Td ^ ri^'^'' [23]. To 
check point (ii), we note that although the quasicondensate is significantly depopulated for 
temperature T > 0.5 Ti^^\ the mechanical effect on the density profile is only noticeable for 
T > 0.8 Tc^^^ [25]. Coherence properties of the strongly degenerate part of the cloud are still 
described by Eq. Hll(l . and provided the correction to ni(0) are taken into account, we expect 
that the scaling behavior still holds to a good approximation, since the correlation function is 
largely insensitive to the precise functional form of ni(z). 

In conclusion, we have investigated in this paper the crossover from a very elongated, 
3D Bose gas to a ID situation where transverse motion is frozen in the limit of vanishing 
density fiuctuations, i.e. for a gas strongly in the degenerate regime. By relying on a local 
density approximation, we have been able to compute the radially integrated density profile 
for any transverse confinement; we believe these results are simple enough to prove useful 
for the analysis of time of fiight images of very elongated samples, with eint ~ hujx. We 
have applied the method to the problem of phase fiuctuations arising in such geometry at 
finite temperatures, and have found an almost "universal" behavior of quasicondensates in an 
elongated geometry, related to the essentially classical nature of thermal phase fluctuations. 
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Table I - Realistic experimental parameters illustrating the various regimes considered in the paper. 
Notations are explained in the text. 





Atom number 






X 








3D regime: 


10" 


1 kHz 


20 Hz 


90 


120 nK 


1.2 pK 


1 mK 


■^•^Na, optical trap 
















Crossover regime: 


5 X 10" 


760 Hz 


5 Hz 


3 


35 nK 


250 nK 


120 AtK 


*^Rb, magnetic trap [12] 
















ID regime: 


200 


20 kHz 


50 Hz 


0.02 


4 nK 


720 nK 


480 nK 


^■^Rb, 2D optical lattice [3] 


















a ni z [az2/a^] 

Fig. 1 - Accuracy of the local density approximation, (a) Local chemical potential near the 3D- 
ID crossover, as a function of the local ID density n\. The circles show the results of a numerical 
calculation [18], undistinguishable at the scale of the figure from Eq. Q [solid line]. The dotted and 
dashed lines show the 3D and ID Thomas-Fermi limiting cases, (b) Integrated density profiles taking 
the (harmonic) axial trapping potential into account. To highlight the crossover, the parameter x = 1 
has been chosen (corresponding to ~ 1.85fiaj_L). The circles, resulting from a numerical solution of 
the Gross-Pitaevskii equation, are indistinguishable from the LDA result (solid line) at the scale of 
the figure. The dotted and dashed lines give the 3D and ID Thomas-Fermi profiles, extrapolated to 
X = 1 for comparison. 
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Fig. 2 - Spatial correlation function of an elongated, phase-iiuctuating condensate for various con- 
finement regimes. In (a), the spatial correlation function is drawn in the 3D case as a function of 
the reduced distance s/L, for T = (upper curve) and for T = (lower curve). The solid line 
follows from a numerical calculation based on the results in [8], and the dashed line is the LDA. 
Figure (b) shows the corresponding curves for the ID ease [5]. In (c), the correlation function of an 
elongated condensate is plotted for T = 4T^ and different regimes of transverse confinement: x = 100 
(3D case, dotted), x = 1 (intermediate case, solid) and x = 10"^ (ID case, dashed). Despite 4 or- 
ders of magnitude of variations in x, the functional form of the spatial correlation function is almost 
unchanged. 



